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B 

' ^ ' We study transport through double quantum dots coupled to normal and superconduct- 

5-H ' ing leads, where the Andreev reflection plays a key role in determining characteristic 

O ' transport properties. We shall discuss two typical cases, i.e. double dots with serial 

vJ , or parallel geometry. For the parallel geometry, the interference of electrons via multi- 

ple paths is induced, so that the transmission probability has Fano-type dip structures 
which are symmetric with respect to the Fermi energy. We also investigate the Aharonov- 
Bohm(AB) effect for the parallel geometry. In some particular situations, we find that 
the general AB period for double dots. An, is reduced to 2it. 



(N 

^D ' 1. Introduction 

^: 

f^ , Recent advances in nanotechnology have enabled us to realize mesoscopic normal- 

metal/superconductor hybrid systems. In these systems, the Andreev reflec- 
tion plays an important role for quantum transport. In particular, the An- 
dreev reflection for a quantum dot coupled to normal and superconducting leads 

'■^ ' gives rise to chara cteristic tra nsport properties due to the discreteness of en- 

C , ergy levels in a dot! I l ''^ ^ " ' ' ! '^ I Moreover, the interplay between the Andreev re- 

flection and the Kondo effect in quantum dot systems has been investigated 
intensively.l3HEIininiTlTl 

In this worH^, we study transport through double quantum dots (DQD's) cou- 
/\ ' pled to normal and superconducting leads. Here, we concentrate on transport due 

jrt ' to the Andreev reflection, which is referred to the Andreev tunneling in the follow- 

ing. We shall discuss two typical cases: the DQD is connected in series or parallel. 
In particular, we focus on the interference effect, which is caused via multiple paths 
of electron propagation, on the Andreev reflection in the parallel DQD . '^ " ' 
Since the interference effect is sensitive to the magnetic flux, we also investigate the 
influence of the Aharonov-Bohm(AB) effect on the Andreev tunneling. 

In the following, we first give a brief explanation of the model, and then describe 
the results for the differential conductance due to the Andreev tunneling in Sec. |31 
A brief summary is given in the last section. 
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2. Model 

We consider a DQD system coupled to normal and superconducting leads 
(N/DQD/S) shown in Fig. Q] In this figure, two dots are coupled via the inter-dot 
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Figure 1. DQD system coupled to normal(N) and superconducting(S) leads. 4c is the inter-dot 
tunneling, and F ^ •^ (i = 1, 2) represents the resonance width due to the transfer between the dot-i 
and normal(superconducting) lead. Note that two dots are connected in series when F^ 
(the dashed arrows), which corresponds to a = (see text). 
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tunneling tc, and the dot-i is connected to the normal (superconducting) lead via 
the tunneling, which causes the resonance width of F^ ^ {i = 1, 2). Introducing the 
ratio a — T^ /T^ — Ff /Ff in the same notation as in Refs. I17I19L we discuss two 
typical cases a = and a ^ I, namely, two dots are connected in series or parallel. 
We assume that the superconducting lead is well described by the BCS theory with 
a superconducting gap A. In addition, the intra-dot Coulomb interaction is ignored 
for simplicity. Since we are interested in the Andreev tunneling, we concentrate 
on the zero temperature case (T — 0) in the region of small bias voltage V (i.e. 
\V\ < A). We calculate the differential conductance dl/dV {I: current) and the 
density of states (DOS) of the dots by using the Keldysh Green functions. 



3. Numerical Results 

3.1. Andreev tunneling in serial and parallel DQD systems 

We first discuss the Andreev tunneling in the serial and parallel DQD systems. For 
simplicity, we fix the energy level of the dot-i (e^) at the Fermi energy (ei = £2 =0), 
and use the gap A as the unit of energy. 

Figure |21 shows the differential conductance as a function of the bias voltage V 
for (a) serial (a = 0), (b) parallel {a — 0.7) DQD system. In both DQD systems, 
the conductance has four peaks in its voltage dependence, which are symmetric 
with respect to the Fermi energy. We shall discuss characteristic aspects of these 
peaks in terms of the DOS of the dots. Note first that the interdot coupling tc 
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Figure 2. The conductance as a function of the bias voltage V for the N/DQD/S system, (a) 
serial DQD (a = 0), (b) parallel DQD (o = 0.7). We set tc = 0.2, F^ = 0.01, Ff = 0.1 (A=l). 



forms the bonding and anti-bonding states for electrons in dots, whereas those for 
holes are obtained by inverting the DOS profile with respect to the Fermi energy. 
The Andreev reflection at the DQD/S-lead interface mixes these states for electrons 
and holes, giving rise to the four Andreev bound states in the dots. Therefore, as 
mentioned above, the DOS of the dots has the four peaks, where the width of these 
peaks is determined by the resonance width Tf' . In the serial DQD system, these 
peaks are indeed observed in the profile of the conductance as shown in Fig. |2Ia). 
On the other hand, in the parallel DQD, Fano-type dip structures {V ~ ±0.2) 
appear in addition to the four peaks, as shown in Fig. IHb). Here, it is instructive to 
recall the interference effect in the parallel DQD system coupled to two normal leads 
(N/parallel-DQD/N). In this system, the DOS for electrons in the dots consists of 
sharp and broad resonance peaks, as shown by the thick line in Fig. I^b). When 
an electron transports via these resonances, its transmission probability acquires a 
Fano-t ype dip structure around the position where the sharp peak in the DOS is 
We note that the DOS for holes is obtained by inverting that for 
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Figure 3. (a) DOS of the dot for the N/parallel-DQD/S system. The solid (dashed) line is for 
the dot-l(2). The parameters are the same as in Fig. |3b). (b) Sketch of DOS for the N/parallel- 
DQD/N system. The thick (thin) line is for electrons (holes), e = corresponds to the Fermi 
energy. 



electrons with respect to the Fermi energy (thin line in Fig. Etb)). Coming back 
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to the N/parallel-DQD/S system, we now see that mixing of the electron and hole 
states induced by the Andreev reflection gives rise to four Andreev bound states 
in the dots, similarly to the serial case. In contrast to the latter case, however, 
the resonances in the parallel case have two different widths: two sharp peaks 
and two broad peaks shown in Fig. El^a). The interference between the distinct 
transport channels via these resonances gives rise to the Fano-type dip structures 
in the conductance, which are symmetric with respect to the Fermi energy (Fig. 
Hb)). 

Here, we make a brief comment on the special case a = 1, i.e. the symmetric 
couplings with leads. In this case, the sharp peaks in Fig. (SJa) become of delta- 
function type, which means that the corresponding local states in the dots are 
completely decoupled from the leads. Therefore, the remaining states with the 
broad resonances in Fig. ISfa) only contribute to the Andreev reflection, so that the 
transport shows similar behavior to the case of a single dot coupled to normal and 
superconducting leads. 



3.2. AB effect in N/DQD/S system 

We next discuss how the Andreev tunneling changes its character when the magnetic 
flux is added in the parallel DQD system. Here, we assume that the magnetic flux 
equally pierces the two subrings formed by the interdot tunneling ic, so that the 
effect of the magnetic flux is symmetrically incorporated in the tunneling between 
the dot and the lead. 

Before considering the N/DQD/S system, we briefly mention the AB effect in 
the N/DQD/N system. As noted in the recent literature,'2ii'2J-22^ ^he interdot 
coupling between the dots forms two-subring structure, so that the AB period in 
the N/DQD/N system becomes 47r instead of the normal period of 27r. Similarly, 
the AB period in the N/DQD/S system is expected to be 47r in general, as pointed 
out by Zhang et al.\^ We find, however, that the AB period is reduced to 27r in 
some particular situations, as will be explicitly shown below. 
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Figure 4. The conductance as a function of V for several choices of the magnetic flux </>. The 
other parameters are the same as in Fig. |3b). 
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Figure 21 shows the conductance for several values of (f>, where (f) represents the 
normalized value of the magnetic flux $ as = 27r<i>/$o (*&o = h/e). The other 
parameters are the same as in Fig. |2Ib). It should be noted that the conductance 
is symmetric with respect to the Fermi energy for any value of (j) and has the AB 
period of 27r (not An). More precisely, the conductance at the magnetic flux 0, 
G{4)), satisfies 0(0) = G{2tt — (f). Including the situation for Fig. 21 we obtain the 
general condition that reduces the AB period to 27r, 

(i) £1=62 = 0, (0<a< 1) 
(ii) El = —£2(7^ 0) and a = 1. 

The condition for the energy levels of the dots, £1 = —£2, means that the bonding 
state and the anti-bonding state in dots are symmetric with respect to the Fermi 
energy [e = ±. e'^t^s + <c )■ Then, the DOS of the dot- 1(2) at the magnetic flux (j), 
Pi(2)(0) has the AB period of 27r {pii^2){4>) = Pi(2)(0 + 27r)) for the case (i). On the 
other hand, Pi{2){4>) for the case (ii) has the periodicity of Pi(2){4') = P2{i){4' + 27r), 
namely, the DOS of the dot-l(2) in the absence of the magnetic flux is the same 
as that of the dot-2(l) for = 27r. Since the DOS with this periodicity directly 
determines electron transport, the conductance changes in the AB period of 27r. 

In this connection, we stress the diffe rence in the AB period between our result 
and the related work by Zhang et a/M^, who treated the same DQD system. We 
have found here that the AB period is reduced to 2tt in the case (ii) , although they 
claimed that it still remains 47r in the same condition. 

4. Summary 

We have studied transport properties through the DQD coupled to normal and 
superconducting leads. It has been discussed that the four-peak structure in the 
DOS is formed by the bonding and antibonding states in the dots coupled to the 
superconducting lead. This structure in the DOS indeed determines the charac- 
teristic properties in the differential conductance, in accordance with the results of 
Zhang et al.W^lw particular, in the parallel DQD system, the interference effect due 
to the multiple paths gives rise to the Fano-type dip structures in the differential 
conductance, which are symmetric with respect to the Fermi energy. We have also 
found the interesting fact that the AB period is reduced to 27r in some particular 
situations, which is contrasted to the AB period 47r expected generally for DQD 
systems. 
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